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^" ■ Abstract 

>. 
(~ — , Masses of singly and doubly charmed baryons are calculated in quenched 

lattice QCD using an improved action of the D234 type on an anisotropic 

j — . ■ lattice. The mass differences between spin 3/2 and spin 1/2 baryon states 

are calculated and compared to mass differences between vector and pseu- 

f— n . doscalar mesons. The suppression of spin splittings in mesons containing 

"pq | heavy quarks, characteristic of quenched QCD simulations, is not observed 

Qh| in the baryon sector. The mass dependence of color hyperfine effects is 

O ,' discussed within the context of the quark model and heavy quark effective 

theory. 
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I. INTRODUCTION 

Quenched lattice quantum chromodynamics (QCD) does reasonably well at describing 
hadronic phenomena. However there are a few well established instances where quenched 
approximation clearly fails to reproduce experimental values. One well studied example is 
the persistent underestimate of mass differences between vector and pseudoscalar mesons 
containing heavy quarks [|| . While the situation is still unsettled, the first attempts with 
unquenched simulations do not yet resolve this discrepancy 0-§| so, clearly, a better 
understanding of color hyperfine effects in lattice QCD would be helpful. One way to 
proceed is to explore the situation in heavy baryons. Unfortunately, this is not so easy 
since spin splittings in the baryon sector are smaller than in mesons while lattice QCD 
correlators for baryons (especially spin 3/2 states) are noisier than those of mesons. An 
early work || suggested that spin splittings in heavy baryons are indeed very small, but 
later studies for charmed baryons || and for bottom baryons @ obtained results more 
in line with phenomenological expectations. 

In this paper we present further results for charmed baryons using the approach of 
H . Calculations were carried out for lattices with different values of the lattice spacing. 
As well, the physical size of the lattice was increased compared to that used in ||. 

Section II presents some details of the lattice simulation. The calculations are done 
with an improved action of the D234 type || on an anisotropic lattice. This action 
has been used and described previously so the detailed expressions are relegated to an 
Appendix. The results are given in Section III along with some discussion about the 
fitting procedure and about the estimate of systematic errors. No attempt was made in 
the present work to do a continuum extrapolation. At each lattice spacing, a number of 
systematic effects have been identified and estimates were made for uncertainties induced 
in masses and mass differences. These include the overall scale uncertainty, the choice of 
time window for fitting correlation functions, ambiguity in fixing the strange and charm 
quark mass and in the extrapolation to physical up and down quark masses. 

In section IV we try to compare spin splittings in baryons with those in mesons. It 
turns out that, for a subset of mesonic and baryonic states, there is a very simple mass de- 
pendence of color hyperfine effects. In fact over the whole mass range where experimental 
masses are available the ratio of meson to baryon mass difference is remarkably constant. 
This fact is used as a benchmark against which to view our lattice QCD results. It is 
suggested that, whereas quenched lattice QCD underestimates mesonic spin splittings, 
in baryons the situation is qualitatively different with no suppression of spin splittings 
being observed. Rather there is a tendency for spin splittings to be overestimated. 

II. METHOD 

The calculation is done on an anisotropic lattice using gauge field and quark actions 
that are improved || by removing tree level errors up to 0(a 2 ) and introducing renor- 
malization factors for the gauge links to remove the dominant errors due to tadpole 



fluctuations. The tadpole factors are estimated using the expectation value of the gauge 
field link in Landau gauge. These actions have been used previously for heavy-light 
mesons || and in a preliminary study of charmed baryons ||. The expressions for the 
actions are given in the appendix. 

Hadron masses are calculated from zero-momentum correlation functions in the usual 
way. This choice of interpolating operators for the hadrons is not unique. For S-like 
baryons containing quarks with two different flavors (denoted by Q and q) a common 
choice |1(J is 



e abc [qIC l5 Q b ]q c , (2.1) 

where a,b,c are color indices and Dirac indices have been suppressed. Taking q = u and 
Q = d, this operator will give the usual interpolating operator for the proton and, as 
discussed in Ref. 0, it is advantageous to use this operator for all flavor combinations 
as this gives a unified analysis over the whole mass range. For the Eg which has light 
quarks of different flavor coupled to spin 1, the operator 

-L {e abc [q' a T C l5 Q b }q c + e abc [q T a C l5 Q b ]q' c } (2.2) 

symmetrized in the light quarks is used. The A-like baryons containing a heavy quark 
can be interpolated using an operator of the form 

e abc [qIC l5 q' b ]Q c . (2.3) 

A more symmetrical choice would be the octet lambda 

-±=e abc {2[qJC l5 q' b }Qc + [qIC l5 Q b ]q' c - [q' a T C l5 Q b }q c } (2.4) 

which is degenerate with the £ in the SU(3) flavor limit. This choice is consistent with 
the idea of using operators that allow for a unified calculation over all masses. In fact, 
it was found that the operators ( |2.3| ) and ( |2.4| ) give masses which are consistent within 
statistical errors so results using (|2.4j ) are reported in this paper. 

For spin 3/2 states we use the simple operator (appropriate for the heavy quark limit) 

e abc [qIC lliq ' b }Q c (2.5) 



for all states since there are indications |T| that the more symmetrical form used in || 
leads to essentially identical results. 



The spin 3/2 field ( |2.5|) propagates both spin 1/2 and spin 3/2 states [12 . At zero 



momentum the correlation function with spatial Lorentz indices has the general form || 
Cij{t) = (<% - \^l 3 )C- i/2 {t) + ± iaj C l/2 (t), (2.6) 



where the subscripts 3/2 and 1/2 denoted the spin projections. By choosing different 
Lorentz components the quantity Cs/2(t) is extracted and used to determine the mass of 
the spin 3/2 states. 

Hadron correlators were calculated using interpolating operators in local form at both 
source and sink and also applying a gauge invariant smearing to the quark propagators at 
the sink. The gauge invariant smearing function, Eq. (13) of Ref. ||13|| , was used. Hadron 
masses were obtained by a simultaneous fit to local and sink-smeared correlators taking 
into account correlations between different time slices with the inverse of the covariance 
matrix obtained using singular value decomposition. Two exponential functions are used 
in the fit to the local correlator. The sink-smeared correlation function is fit by a single 
exponential using a mass parameter constrained to be the same as the ground state mass 
of the fit to the local correlator. The time window for the fit was obtained by fixing the 
ending time to be sufficiently large so that the fits were not sensitive to its value and 
by varying the starting time so that a minimum of x 2 (typically around 1/d.o.f.) was 
achieved. Fits with starting times ±1 time step away from the time giving the minimum 
X 2 were used to estimate the systematic uncertainty associated with fixing the fitting 
time window. 

III. RESULTS 

Calculations were carried out for three different sets of quenched gauge configurations 
on anisotropic lattices with a bare aspect ratio a s /a t = 2 and spatial lattice spacing vary- 
ing from about 0.22fm to 0.15fm. Gauge fields were constructed using pseudo-heatbath 
Monte Carlo with 400 {(3 = 2.1) to 800 (/? = 2.5) sweeps between saved configurations. 
Fixed time boundaries were used in calculating the quark propagators. The parameters 
for the lattices are given in Table |I[ 

Six values of the hopping parameter were used and are listed in Table |TJ The smallest 
values were chosen in the region of the charm quark and the four larger values were used 
for extrapolating to the light (u,d) quark mass. The temporal lattice spacing was fixed by 
calculating the p-meson mass. As well known, quenched lattice QCD does not give a mass 
spectrum in perfect agreement with experiment so alternate ways to fix the scale would 
lead to different values. This is an intrinsic limitation of the quenched approximation. 
The hopping parameters corresponding to the strange and charm quark masses were fixed 
from the 0-meson and D-meson respectively. These are the values given in Table [ll[ A 
systematic uncertainty in these quark mass determinations was estimated by also using 
the kaon and J/ip masses. 

The masses of hadrons containing up and down quarks have to be extrapolated into 
the light quark mass region. This is done by extrapolating the masses calculated at the 
four largest hopping parameters as a function of pion mass using the function cq + c^vp?^ + 
czm%. In some cases the mass to be extrapolated appears to be described very well by a 
function without an m^ term. In such cases, where the coefficient C3 comes out to be not 



statistically significant, the difference between aquadratic extrapolation and a quadratic 
plus cubic extrapolation is included as a contribution to the systematic error. 

Some representative results of the mass determination are shown in terms of so-called 
effective mass plots for singly heavy baryons in Figs. 1 and 2. The effective mass 
M(t) is ln(g(t)/g(t + 1)) where g(t) is the zero-momentum time correlation function for 
the hadron. The solid line in the figures shows the ground state mass obtained by a 
simultaneous fit to the local and smeared correlation functions plotted over the time 
window of the fit, and the dashed lines indicate the bootstrap error. Note that a two 
exponential form was used to fit the local correlators. Corresponding results for doubly 
heavy baryons are given in Figs. 3 and 4. 

The results for singly charmed baryons are given in Table |T|. The first error is the 
statistical error which is calculated using a bootstrap analysis employing a bootstrap 
sample size equal to the configuration sample size. The second error is the combined 
systematic error. This includes the overall scale uncertainty, the uncertainty due to 
choice of the correlation function fitting range as well as strange and charm quark mass 
uncertainties and light (u,d) quark mass extrapolation ambiguity, where applicable. For 
comparison, the experimental values are also shown, where they are known. 

The masses of doubly charmed baryons were also calculated and are given in Table 
P"V| . No doubly heavy baryons have been observed yet in experiments but one may hope 
that this situation will change in the not too distant future (see Ref. [0 for a review of 



the possibilities for experimental observation). Note also that the spin splittings for the 
doubly charmed baryons are as large and perhaps even larger than in the singly charmed 
sector. 

IV. DISCUSSION 

In order to see how our results for charmed baryons and particularly for the spin 
splittings fit into the overall scheme of hadron physics, it is useful to start with a slight 
digression and consider the spin splittings of pseudoscalar and vector mesons. It is well 
known that for vector (V) and pseudoscalar(P) meson pairs of the form Qq, where q is 
up or down and Q is any light or heavy flavor, the squared mass difference M v — Mp is 
approximately constant for all flavors Q. This mass relation was first derived within the 
framework of string-like models for hadrons imposing the constraints of chiral symmetry 
i~7HT9|. It can also be understood in the quark model with a linear confining potential 



20fl and, for heavy-light mesons, from heavy quark effective theory (HQET) [^TJ. For 
our purposes it is useful to factorize the squared mass difference and make a plot of 
the spin splitting AM mes = My — Mp versus the inverse of the average meson mass 
M ave = {My + Mp)/2. This is given in Fig. 5 where the line is the best linear fit to 
the experimental points shown by triangles. The mesons pairs included in the plot are 
(p, 7i"), (K*, K), (D*, D) and (B*, B). The remaining points in Fig. 5 are the results from 
quenched lattice QCD calculations. The squares are results of the present simulation 
and the circles are representative results collected from a variety of published papers 



|7|J9| J2^J23|1 . The underestimate of color hyperfine effects by quenched lattice QCD can be 
seen clearly 

We would like to have a similar global view of spin splittings in the baryon sector. 
In fact the string models that predict the meson squared mass difference relation give 
a similar relation for baryons [TH]. However, if one plots the baryon mass difference 



AM{, ar = M3/2 — M1/2 versus the inverse of the average baryon mass (M 3 / 2 + M 1 / 2 )/2 
for the well measured baryon pairs (A, N), (E*,E) and (E*,E C ), it is found empirically 
that the relation is not linear. On the other hand, if one uses the average meson mass 
(associating mesons and baryons with the same flavor content) an almost exact linear 
relation is found. This is shown in Fig. 6 where the triangles are experimental data and 
the line is the best linear fit. The results of the present simulation as well as some results 
of published quenched lattice calculations [p|P|] are also shown in Fig. 6. Unfortunately, 



masses in the baryon sector are not determined with the same precision as for mesons 
but a clear qualitative difference from Fig. 5 can be seen. No suppression of the spin 
splitting is evident. Although not completely conclusive, there may be a tendency for 
baryon spin splittings to be overestimated. 

The linearity of the experimental results for mass differences plotted versus a common 
variable in Figs. 5 and 6 implies that the ratio of mass differences should be constant. 
This is shown in Fig. 7 for the experimental data where the ratios of the meson to baryon 
mass differences (p, vr)/(A, iV), (K*,K)/(E*,E) and (£>*, D)/(E$, E c ) are plotted versus 
the averaged meson mass. The ratio is constant to within about 2% and the solid line is 
the average value 2.13. This remarkable result was anticipated by Lipkin ]25| (see also 
Lipkin and O'Donnell PEfl ) within the framework of a quark model although it required 
a number of assumptions. For mesons and baryons with a single heavy quark, HQET 
also implies a constancy in the meson to baryon spin splitting but the effective theory 
can not predict the value of the ratio. 

One would like to see what lattice QCD predicts for the meson to baryon mass 
difference ratio. Unfortunately the errors associated with extrapolation preclude a very 
precise determination. For this reason we choose to present unextrapolated results, fixing 
the "light" quark at a kappa value near that of the strange quark. The results for 
Qqq — baryons (q-mass fixed, Q-mass variable) are shown in Fig. 8. The calculated ratio 
is quite constant especially for the simulations at higher (3 where the results are more 
precise. The fact that the ratio falls below the empirical value, shown by the solid line in 
Fig. 8, is an indication that the suppression of spin splittings is not present in baryons 
in the same way as it is in mesons. 

For doubly heavy baryons one can make definite predictions for the relation between 
mesonic and baryonic spin splittings. Consider the simplest possible quark model in 
which only quark masses and the color hyperfine interaction term are taken into account 
2T| , |2"8"|| . The operators for meson and baryon mass are 



M mes = m 1 + m 2 + a(al-a^)/m 1 m 2 , (4.1) 

and 



M bar = m 1 + m 2 + m 3 + - ^ a^(ot ■ Oj )/m i m j . (4.2) 



Note the factor of 1/2 in the last term of ( |4.2|) . This reflects the reduction in the strength 
of the gluon exchange between quarks in a color anti-triplet state relative to that between 
a quark and antiquark in a color singlet. Evaluating the above expressions for mesons 
Qq and H baryons QQq one finds 

3r 

M(P)=m q + m Q , (4.3) 

m q m Q 

M{V)=m q + m Q + —^—, (4.4) 

m q m Q 



M(E QQ ) =m g + 2m Q + 2\-^- —S3- , (4.5) 

{ 4m z Q m q m Q J 

M(E* 00 ) =m q + 2m Q + 2 I ^- + ° qQ 1 , (4.6) 

where the coefficients c and d depend on the probability of finding the interacting quarks 
at zero separation. The relation between the spin splittings is 

AM bar = ^AM mes . (4.7) 

4 c 

Now, if it is assumed that the two heavy quarks act as a single heavy compact antitriplet 
color source, it is reasonable to expect that d qQ ~ c so that 

AM bar w -AM mes . (4.8) 

This diquark picture for doubly heavy baryons can also be analyzed in a heavy quark 
effective theory. This has been done in an elegant formulation utilizing a superflavor 
symmetry to relate hadrons containing a heavy vector diquark to those with a heavy 



spin-1/2 antiquark [29 



Although no evidence yet exists that for doubly heavy baryons the ratio of meson to 
baryon spin splittings is a constant function of quark mass our lattice calculations suggest 
that it might be so. Fig. 9 shows the ratio of unextrapolated results for QQq — baryons 
(q-mass fixed, Q-mass variable). The long dashed line in the figure is the value of the 
ratio obtained experimentally for the strange quark, that is the K* — K mass difference 
divided by the S* — S mass difference. The short dashed line is the prediction Eq. ( |4.8| ) . 

Eq. ( f4.7|) has been derived explicitly for the doubly heavy Hqq, Sqq system but 
the same expression holds for the singly heavy baryons. In this case there is no simple 
relation between c'q and c but, since the gluon interaction between quarks is weaker than 
between a quark and antiquark, one might expect that d qQ < c. This implies that spin 



splittings in singly heavy baryons are smaller than in doubly heavy baryons containing 
the same quark flavors. This expectation is borne out empirically for the strange E's 
and Efs but given the statistical and systematic errors there is no clear evidence from 
our lattice simulations for charmed baryons. 

Finally, we note that if it is assumed that the ratio AM mes /AM bar is constant as a 
function of quark mass for doubly heavy baryons just as it is for singly heavy baryons, 
then using the experimental values of the S and 5* masses along with the known meson 
masses yields the phenomenological predictions 

~* cc - ~ cc = 76.6MeV, 

and 

E* bb - ~ bb = 24.5MeV. 

V. SUMMARY 

Masses of charmed baryons were calculated in quenched lattice QCD using improved 
gluon and quark actions on an anisotropic lattice. The actions were improved to re- 
move tree level errors up to 0(a 2 ) and tadpole factors, estimated by using gauge field 
links in Landau gauge, were introduced to remove the dominant errors due to tadpole 
fluctuations. 

Calculations were done at three different values of the gauge coupling constant with 
the spatial lattice spacing varying from 0.22fm to 0.15fm. The results at the two largest 
values of the gauge coupling are compatible with each other and are in fair agreement 
with experimental data where it is available. Masses and mass differences at the smallest 
value of the gauge coupling (largest lattice spacing) tend to be smaller, perhaps reflecting 
larger finite lattice spacing errors. 

The main focus of our work was the mass differences due to spin dependent inter- 
actions. It is well established that in quenched lattice QCD the mass differences be- 
tween vector and pseudoscalar mesons are underestimated for mesons containing heavy 
quarks. The present simulations show no comparable suppression of splittings between 
spin 1/2 and spin 3/2 baryons. Our results and results taken from the literature indicate 
a tendency for baryon spin splittings to be overestimated but, this is not established 
definitively due to relatively large errors associated with baryon mass determinations. 

To get an overall view of the spin splittings it was useful to consider the scaling of 
mass differences with the average meson mass as one changes quark flavors. For mesons 
this relation [17[ predates QCD. We have found it useful to extend the same scaling (with 



average meson mass) to mass differences between baryons. A relation between meson and 



baryon spin splittings, which is implied by this scaling, was anticipated by Lipkin [2S| 
from a quark model analysis. From the point of view of QCD, a mass independent ratio 
of meson to baryon spin splittings can be derived for heavy flavored mesons and baryons 

8 



using heavy quark effective theory, but how to extend this result to light flavor hadrons 
is not clear. Nonetheless a constant meson to baryon mass difference ratio is satisfied 
remarkably well by experimental data for quark flavors from light to charm. Obviously, 
a definitive experimental determination of the ££ mass to extend this analysis to the 
b-quark region would be extremely useful. 

Our lattice QCD results are compatible with a mass independent meson to baryon 
spin splitting ratio for both singly and doubly heavy baryons. However, the more precise 
values obtained at our larger values of the gauge coupling indicate a value for the ratio 
for singly heavy baryons smaller than the empirical value. This can be interpreted as 
another indication that the suppression of spin splittings found for mesons does not occur 
for baryons. For doubly heavy baryons our preferred values for the meson to baryon ratio 
lie below the value obtained by using the masses of the doubly strange hyperons S and 
H*. Clearly, the experimental observation of doubly heavy baryons and the systematic 
investigation of their spin splittings would be very interesting indeed. 
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APPENDIX A: THE LATTICE ACTION 

The lattice action has two terms: gauge action and quark action. The entire action 
is classically and tadpole-improved with the tadpole factors, u s and u t , defined as the 
mean links in Landau gauge in a spatial and temporal direction, respectively. 

The gauge field action is 



S G (U) 



5(3 

T 



-Ls^-iReT^-^S^-iReW,, 






(Al) 



where £ = a s /a t is the aspect ratio and (3 is the lattice gauge field coupling constant. 
The subscripts "ps" and "rs" denote spatial plaquettes and spatial planar 1x2 rectan- 
gles respectively. Plaquettes in the temporal-spatial planes are denoted by "pt", while 
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rectangles with the long side in a spatial(temporal) direction are labeled by "rst"("rts"). 
The leading classical errors of this action are quartic in lattice spacing. 

An action of the D234 type || is used for the quarks with coefficients set to their 
tadpole- improved classical values. Its leading classical errors are cubic in lattice spacing: 



S F (q,q;U) 



4k 



3 
+ 



E 

x,i 

4k 



1 



E 



I 

Ut 



D u (x) - 



D u (x) 



1 



1 



D 2i (x) 



to?" 



D 2t (x) 



2k 



u-2 E ^{x^ijFijix)^) 



Me i 

2k 
3u 2 a uU 



Kj 



^2ip(x)a 0i F 0i (x)ij(x) 



E^OM 



X 



(A2) 



where k denotes the hopping parameter and 


D u (x)=$(x)(l 


-^Ti) c/ i(a ; )V'(^ + «) 


+ 4>(x + 


i)(l + £ 7i )a?(sMx), 


D lt (x)=$(x){l 


— l±)U A (x)ip{x + i) 


+ 4>(x + 


?)(1 + 74)^4(3:)^), 


D 2l (x) = 4>(x)(l 


-^7 i )C/ i (x)C/ i (xH-z)V(^ + 2l) 


+ 4){x + 2z)(l + in)U}{x + i)Ul(x)iJ)(x), 


D 2t {x)=4>{x){l 


- 7 4 )f/4(x)f/ 4 (x + £)V>(z + 2£) 


+ 4>{x + 


2i)(l + 7 4 )C/ 4 t (^ + t)t/ 4 f (x)^(a;), 


gF, v {x) = i (^(x) - ^( x )) _ il m (THV(x)) 



"/i!/ 



H,,,, = — [f/ M (a;)^(x + //)C/t(z + J>)E#(x) 

+t/„(a?)L^(x - /t + v)Ul{x - fi)U^(x - ft) 
+Uj l (x — jj)Ul(x — fi — u)Un(x — fi — v)U v (x — 0) 
+Ul(x-v)U li {x-i>)U u {x + fi-v)Ul{x) . 



(A3) 
(A4) 
(A5) 

(A6) 

(A7) 



(A8) 
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TABLES 







TABLE I. Parameters for lattices used in 


;his work. 




p 


size 


configurations 


ai\GeV) 


u s 


Uf 


2.1 


12 3 x 32 


720 


1.803(42) 


0.7858 


0.9472 


2.3 


14 3 x 38 


442 


2.210(72) 


0.8040 


0.9525 


2.5 


18 3 x 46 


325 


2.625(67) 


0.8185 


0.9564 



TABLE II. Hopping parameter values. The quantities k s and k c are the hopping parameter 
values associated with strange and charm quarks respectively. 







K, S 



K s (4>) 



K C (D) 



2.1 
2.3 
2.5 



0.175,0.176,0.229,0.233,0.237,0.240 
0.183.0.189,0.229,0.233,0.237,0.240 
0.193,0.197,0.230,0.234,0.238,0.240 



0.2338 


0.1739 


0.2371 


0.1875 


0.2382 


0.1964 



TABLE III. Masses of singly charmed baryons. 
are in MeV. The first error is the statistical error 



Masses are given in GeV, mass differences 
and the second is the combined systematic 



error. The experimental values are taken from [14 except 


for E/ c which is from 


Hi- 




[3 = 2.1 


/3 = 2.3 
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69(7) (|) 
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TABLE IV. Masses of doubly charmed baryons. Masses are given in GeV, mass differences 
are in MeV. The first error is the statistical error and the second is the combined systematic 
error. 

(3 = 2.1 /? = 2.3 /3 = 2.5 



cc 



a* 



n* cc -n c 



3.608(15) (i) 


3.595(12)(2i) 


3.605(12)(ll 


3.666(18)$) 


3.678(15)(i|) 


3.685(14) (}? 


3.747(9)$) 


3.727(9) Q 


3-733(9) (J 8 ) 


3.804(13) (11) 


3.800(ll)(i°) 


3.801(9) (f 4 ) 


58(14)(Jg) 


83(8) (/ ) 


80(10)(?) 


57(8) (10) 


72(5) (|) 


68(5) (1) 
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Figure Captions 

Fig. 1 Effective mass M(t) versus t for spin 1/2 baryons Qqq with (a) n q = 0.233, kq = 
0.183 (b) K q = 0.233, kq = 0.189 (c) K q = 0.237, k q = 0.183 (d) K q = 0.237, k q = 0.189. 
The solid line indicates the ground state mass and the dashed lines the statistical error. 

Fig. 2 Effective mass M(t) versus t for spin 3/2 baryons Qqq with (a) K q = 0.233, kq = 
0.183 (b) K q = 0.233, kq = 0.189 (c) n q = 0.237, k q = 0.183 (d) n q = 0.237, k q = 0.189. 
The solid line indicates the ground state mass and the dashed lines the statistical error. 

Fig. 3 Effective mass M(t) versus t for spin 1/2 baryons QQq with (a) K q = 0.233, k,q = 
0.183 (b) K q = 0.233, kq = 0.189 (c) n q = 0.237, k q = 0.183 (d) n q = 0.237, k q = 0.189. 
The solid line indicates the ground state mass and the dashed lines the statistical error. 

Fig. 4 Effective mass M(t) versus t for spin 3/2 baryons QQq with (a) K q = 0.233, kq = 
0.183 (b) K q = 0.233, kq = 0.189 (c) n q = 0.237, k q = 0.183 (d) n q = 0.237, k q = 0.189. 
The solid line indicates the ground state mass and the dashed lines the statistical error. 

Fig. 5 The mass difference between vector and pseudoscalar meson pairs plotted versus 
the inverse of the average mass. Triangles are experimental values, circles are results of 
quenched lattice calculations taken from the literature |7],|9], |22]J23| and squares the results 
from the present work. The line is the best linear fit to the experimental points. 

Fig. 6 The mass difference between spin 1/2 and spin 3/2 baryon pairs plotted versus 
the inverse of the average mass of vector and pseudoscalar mesons with the same flavor. 
Triangles are experimental values, circles are results of quenched lattice calculations taken 
from the literature @,0] and squares the results from the present work. The line is the 
best linear fit to the experimental points. 

Fig. 7 The ratio of the mass difference between vector and pseudoscalar meson to the 
mass difference between spin 1/2 and spin 3/2 baryon for hadrons with the same flavor 
content plotted versus average meson mass. The triangles are experimental data for 
(p,ir)/(A,N), (K*,K)/(E*,E) and (D*,D)/(S*,E C ). The line is the average value. 

Fig. 8 Lattice simulation results for the ratio of the mass difference between vector and 
pseudoscalar meson to the mass difference between singly heavy spin 1/2 and spin 3/2 
baryon at (3 = 2.1 (triangles), (3 = 2.3 (circles) and (3 = 2.5 (squares). The line is the 
average experimental value from Fig. 7. 

Fig. 9 Lattice simulation results for the ratio of the mass difference between vector and 
pseudoscalar meson to the mass difference between doubly heavy spin 1/2 and spin 3/2 
baryon at (3 — 2.1 (triangles), (3 — 2.3 (circles) and (3 — 2.5 (squares). The long dashed 
line is the experimental value for (K*, K)/(S*, E). The short dashed line is the prediction 
from Eq. (O). 
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